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The current-driven motion of a domain wall (DW) in a ferromagnet with many Bloch lines (BLs) via the 
spin transfer torque is studied theoretically. It is found that the motion of BLs changes the current-velocity 
{j-v) characteristic dramatically. Especially, the critical current density to overcome the pinning force is 
reduced by the factor of the Gilbert damping coefficient a even compared with that of a skyrmion. This 
is in sharp contrast to the case of magnetic field driven motion, where the existence of BLs reduces the 
mobility of the DW. 


Domain walls (DWs) and bubbles^’are the spin tex¬ 
tures in ferromagnets which have been studied inten¬ 
sively over decades from the viewpoints of both funda¬ 
mental physics and applications. The memory functions 
of these objects are one of the main focus during 70’s, but 
their manipulation in terms of the magnetic field faced 
the difficulty associated with the pinning which hinders 
their motion. The new aspect introduced recently is the 
current-driven motion of the spin textures.^’The flow 
of the conduction electron spins, which follow the direc¬ 
tion of the background localized spin moments, moves 
the spin texture due to the conservation of the angu¬ 
lar momentum. This effect, so called the spin transfer 
torque, is shown to be effective to manipulate the DWs 
and bubbles compared with the magnetic field. Magnetic 
skyrmion^’®) is especially an interesting object, which is 
a swirling spin texture acting as an emergent particle 
protected by the topological invariant, i.e., the skyrmion 
number iVsk, defined by 
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with n(r) being the unit vector representing the direc¬ 
tion of the spin as a function of the two-dimensional spa¬ 
tial coordinates r. This is the integral of the solid angle 
subtended by n, and counts how many times the unit 
sphere is wrapped. The solid angle and skyrmion number 
A^sk also play essential role when one derives the equation 
of motion for the center of mass of the spin texture, i.e., 
the gyro-motion is induced by fVsk in the Thiele equation, 
where the rigid body motion is assumed.^’ 

Beyond the Thiele equation,one can derive the equa¬ 
tion of motion of a DW in terms of two variables, i.e., 
the wall-normal displacement q{t, Q, rj) and the wall- 
magnetization orientation angle (see Fig. 1) 



Fig. 1. Schematic magnetization distribution of DW with many 
Bloch lines. 


where and t] are general coordinates specifying the 
point on the DW:®^ 
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Here, ' means the time-derivative. || and T indicate 
the components parallel and perpendicular to the DW 
respectively. M is the magnetization, 7 is the gyro- 
magnetic ratio, and a, A are the energy per area and 
thickness of the DW. u® is the velocity of the conduction 
electrons, which produces the spin transfer torque, a is 
the Gilbert damping constant, and (3 represents the non- 
adiabatic effect. These equations indicate that q and ip 
are canonical conjugate to each other. This is understood 
by the fact that the generator of the spin rotation nor¬ 
mal to the DW, which is proportional to sin ip in Fig. 1, 
drives the shift of q. (Note that ip is measured from the 
fixed direction in the laboratory coordinates.) 

In order to reduce the magnetostatic energy, the spins 
in the DW tend to align parallel to the DW, i.e., Bloch 
wall. When the DW is straight, this structure is coplanar 
and has no solid angle. From the viewpoint of eqs. (2) 
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and (3), the angle ip is fixed around the minimum, and 
slightly canted when the motion of q occurs, i.e., ip = 0. 
However, it often happens that the Bloch lines (BLs) 
are introduced into the DW as shown schematically in 
Fig. 1. The angle ip rotates along the DW and the Neel 
wall is locally introduced. It is noted here that the solid 
angle becomes finite in the presence of the BLs. Also with 
many BLs in the DW, the translation of BLs activates 
the motion of the angle ip, i.e., ip 0, which leads to the 
dramatic change in the dynamics. 

In the following, we focus on the straight DW which 
extends along x-direction and is uniform in 2 -direction. 
Thus, the general coordinates here are {(,r]) = {x,z). 
q{t, X, z) is independent of the coordinates q{t, x, z) = 
q(t), and the functional derivative 6a/6q in eq. (3) be¬ 
comes the partial derivative dajdq. In the absence of 
BLs, we set ip{t,x,z) = ip{t), and 6ajSip in eq. (2) also 
becomes da/dip. Then the equation of motion in the ab¬ 
sence of BL is 


dip 


2 M 7 - 


q — aAip — 


(4) 
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dq 
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ip + aA-^q- PA-^vl , 


(5) 


With many BLs, the sliding motion of Bloch lines along 
DW, which activates ip, does not change the wall energy, 
i.e., 6a/6ip in eq. (2) vanishes.Here, for simplicity, we 
consider the periodic BL array with the uniform twist 
ip{t, X, z) = {x—p{t))/A where A is the distance between 
BLs, which leads to 
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First, let us discuss the magnetic field driven motion 
without current. The effect of the external magnetic field 
T^ext jg (described by the force dajdq = —in 
eqs. (5) and (7). and are set to be zero. In the 
absence of BL, as mentioned above, the phase ip is static 
Ip = 0 with the slight tilt of the spin from the easy-plane, 
and one obtains from eq. (5) 

i = ( 8 ) 

a 

This is a natural result, i.e., the mobility is inversely 
proportional to the Gilbert damping a. ip is determined 
by eq. (4) with this value of the velocity q. 

In the presence of many BLs, eqs. ( 6 ) and (7) give the 
velocities of DW and BL sliding driven by the magnetic 
field as 

q=^^A^H^-\ (9) 

1 + 


I 


p = - 






( 10 ) 


Comparing eqs. ( 8 ) and (9), the mobility of the DW is re¬ 
duced by the factor of o? since a is usually much smaller 
than unity. We also note that the velocity of the BL slid¬ 
ing p is larger than that of the wall q by the factor of 
a. Physically, this means that the effect of the external 
magnetic field mostly contributes to the rapid mo¬ 
tion of the BLs along the DW rather than the motion of 
the DW itself. These results have been already reported 
in refs.^’®’ 

Now let us turn to the motion induced by the current 
u®. In the absence of BL, again we put 1 /) = 0 in eqs. (4) 
and (5). Assuming that there is no pinning force or ex¬ 
ternal magnetic field, i.e., dajdq = 0 , one obtains from 
eq. (5) 
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q = -w_L, 
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( 11 ) 


and eq. (4) determines the equilibrium value of ip. When 
the pinning force dajdq = is finite, there appears a 
threshold current density ('c^)^, which is determined by 
putting g = 0 in eq. (5) as 


^ 2Mj3 ’ 


( 12 ) 


which is inversely proportional to Since eq. (II) is 
independent of Uy, the threshold current density 

is (i)^ = oo. 

In the presence of the many BLs, on the other hand, 
eqs. ( 6 ) and (7) give 
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which is the main result of this paper. From eq. (13), the 
current-velocity characteristic in the absence of both the 
pinning and the external field (da jdq=0) is 
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1 4- ct 
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where the fact a, pi <€. 1 is used in the last step. If we 
neglect the term coming from r)j|, the current-velocity 
relation becomes almost independent of a and /? in 
sharp contrast to eq. (11). This is similar to the univer¬ 
sal current-velocity relation in the case of skyrmion,^^^ 
where the solid angle is finite and also the transverse 
motion to the current occurs. Note that rijl slightly con¬ 
tributes to the motion when a ^ pi, while it does 
not in the absence of BL. Even more dramatic is the 
critical current density in the presence of the pinning 
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Fig. 2. The wall displacement g as a fucntion of t for the DWs 
without BL and with BLs. (a) = 22.0. The inset shows the 

pinning force (b) = 21.0. (c) = 0.0043. (d) v^j_ = 

0.0042. 


{da/dq = When we apply only the current per¬ 

pendicular to the DW, i.e., ujj = 0, putting g = 0 in 
eq. (13) determines the threshold current density as 
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which is much reduced compared with eq. (12) by the 
factor of ^ 1. Note that in eq. (15) is even 

smaller than the case of skyrmion^^) by the factor of 
a. Similarly, the critical current density of the motion 
driven by ujj is given by 
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which can also be smaller than eq. (12). 

Next we look at the numerical solutions of q{t) driven 
by the current perpendicular to the wall under the 
pinning force. We assume the following pinning force: 
(yA/2M)iPP”(g) = u*(g/A)exp [-{q/Af] (see the in¬ 
set of Fig. 2(a)). We employ the unit of A = w* = 
1 and the parameters {a,/3) are fixed at {a, (3) = 
(0.01,0.02). Here, we compare two DWs without BL 
and with BLs. The maximum value of the pinning force 
{jA/2M)F^^^ = 0.429 determines the threshold current 
density (u^)^. as (u^)^. = 21.4 and (u^)^, = 0.00429 in the 
absence of BL and in the presence of many BLs, respec¬ 
tively. In Fig. 2(a), both DWs overcome the pinning at 
the current density v‘]_ = 22.0, although the velocity of 
the DW without BL is suppressed in the pinning poten¬ 
tial. At the current density = 21.0 below the threshold 
value in the absence of BL, the DW without BL is pinned, 
while that with BLs still moves easily (Fig. 2(b)). The 
velocity suppression in the presence of BLs is observed 
at much smaller current density = 0.0043 (Fig. 2(c)), 
and finally it stops at = 0.0042 (Fig. 2(d)). 

All the discussion above relies on the assumption that 


the wall is straight and ip rotates uniformly. When the 
bending of the DW and non-uniform distribution of BLs 
are taken into account, the average velocity and the 
threshold current density take the values between two 
cases without BL and with many BLs. The situation 
changes when the DW forms closed loop, i.e., the do¬ 
main forms a bubble. The bubble with many BLs and 
large \Nsk\ is called hard bubble because the repulsive 
interaction between the BLs makes it hard to collapse 
the bubble.At the beginning of the motion, the BLs 
move along the DW, which results in the tiny critical cur¬ 
rent. In the steady state, however, the BLs accumulate 
in one side of the bubble.Then, the configuration 
of the BLs is static and the Thiele equation is justified 
as long as the force is slowly varying within the size of 
the bubble. The critical current density (u®)^ is given by 
(u®)j, (X FP“/iVsk (iVsk ( 3 > 1): the skyrmion number of 
the hard bubble), and is reduced by the factor of Ngk 
compared with the skyrmion with Ask = ±1. 

In conclusion, we have studied the current-induced 
dynamics of the DW with many BLs. The finite ip in 
the steady motion activated by BLs sliding drastically 
changes the dynamics, which has already been reported 
in the field-driven case. In contrast to the field-driven 
case, where the mobility is suppressed by introducing 
BLs, that in the current-driven motion is not necessarily 
suppressed. Instead, the current-velocity relation shows 
universal behavior independent of the damping strength 
a and non-adiabaticity j3. Furthermore, the threshold 
current density in the presence of impurities is tiny even 
compared with that of skyrmion motion by the factor of 
a. These findings will stimulate the development of the 
racetrack memory based on the DW with many BLs. 

Acknowledgments We thank W. Koshibae for useful discus¬ 
sion. This work is supported by Grant-in-Aids for Scientific Re¬ 
search (S) (No. 24224009) from the Ministry of Education, Cul¬ 
ture, Sports, Science and Technology of Japan. J. I. was supported 
by Grant-in-Aids for JSPS Fellows (No. 2610547). 


1) A. Hubert and R. Schafer, Magnetic Domains: The Analysis 
of Magnetic Microstructures (Springer-Verlag, Berlin, 1998). 

2) A. P. Malozemoff and J.C. Slonczewski, Magnetic Domain 
Walls in Bubble Materials (Academic Press, New York, 1979). 

3) J. C. Slonczewski, J. Magn. Magn. Mater. 159, L1-L7 (1996). 

4) L. Berger, Phys. Rev. B 54, 9353-9358 (1996). 

5) S. Miihlbauer et ah. Science 323, 915 (2009). 

6) X. Z. Yu et ah, Nature 465, 901 (2010). 

7) A. A. Thiele, Phys. Rev. Lett. 30, 230 (1973). 

8) K. Everschor et ah, Phys. Rev. B 86, 054432 (2012). 

9) J. C. Slonczewski, J. Appl. Phys. 45, 2705 (1974). 

10) A. P. Malozemoff and J. C. Slonczewski, Phys. Rev. Lett. 29, 
952 (1972). 

11) G. Tatara et ah, J. Phys. Soc. Japan 75, 64708 (2006). 

12) J. Iwasaki, M. Mochizuki and N. Nagaosa, Nat. Commun. 4, 
1463 (2013). 

13) G. P. Vella-Coleiro, A. Rosencwaig and W. J. Tabor, Phys. 
Rev. Lett. 29, 949 (1972) 

14) A. A. Thiele, F. B. Hagedorn and G. P. Vella-Coleiro, Phys. 


3 
































J. Phys. Soc. Jpn. 


LETTERS 


Rev. B 8, 241 (1973). 

4 



